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Quasiparticle spectrum associated with the nodal structure in d-wave superconductors is of great 
interest. We study theoretically the quasi-particle spectrum in a planar magnetic field, where the 
effect of the magnetic field is treated in terms of the Doppler shift. We obtain the angular dependent 
specific heat in the presence of a planar magnetic field and impurities, both in the superclean 
hmit(.^ <C ^ 1) and in the clean limit(.^^ <C <S 1). Also a similar analysis is used for the 
thermal conductivity tensor within the a-b plane. In particular, in contrast to the earlier works, we 
find a fourfold symmetry term in and kx —Hsm{29) where Ky and it± are the diagonal- and 
the off-diagonal components of the thermal conductivity tensor and 9 is the angle between the heat 
current and the magnetic field. 

PACS numbers: 74.60.Ec, 74.72.-h, 74.25. Fy, 74.25. Bt 



I. INTRODUCTION 

In the last few years c?-wave superconductivity in 
hole-doped high-Tc cuprates is established.BB More re- 
cently, Tsuei and Kirtley and other groupa3 have 
shown by means of the phase sensitive tri-crystal ge- 
ometry and the magnetic penetration depth measure- 
ment that the superconductivity in the electron doped 
high-Tc cuprates NCCO and PCCO is also of d-wave. 
This suggests strongly the universality of d-wave su- 
perconductivity for all high-Tc cuprate superconduc- 
tors. Further, recent measurements of the magnetic 
penetration depth in organic superconductor k-(BEDT- 
TTF)2Cu(N(CN)2)Br indicates this system has also d- 
wave superconductivityjjia though the phase sensitive ex- 
periment is not available for this system. Also rincertain 
is the nodal directions of the order parameterSQ in k- 
(ET)2 salts in spite of th&,tiieoretical suggestion that it 
belongs to da;y-symmetry.El"t3 In this circumstance, the 
experimental probe which indicates the nodal directions 
in d-wave superconductors is of prime importance. 

We recall earlier study of the thermal conductivity ten- 
sor in the optimally doped YBCO in a planar magnetic 
field indicates clearly the nodal directions of d-waV|&.su- 



perconductors consistent wit 



It l^Q 



symmetry.Ert3 



In the meanwhile VolovikEj made the crucial obser- 
vation that the thermodynamics of the vortex state in 
d-wave superconductors is dominated by the extended 
quasi-particle states associated with the nodal structure 
and the effect of the magnetic field is incorporated by the 
Doppler shift in the energy spectpHm.li3 This approajehia 
further developed by Barash et aO and Kiibert et al.ll3't§ 

In particular, those authors consider the supercurrent 
associated with individual vortices explicitly and then 
made spatial average of individual vortex contribution 
over the unit cell of the vortex lattice. Indeed, Kiibert et 
al have succeeded not only to describe the earh2-.specific 
heat data by K. Moler et al and B. Revaz et aMM, but 
also predicted the thermal conductivity in a magnetic 



field parallel to the c-axis in the low temperature limit, 
which is confirmed byj-the thermal conductivity measure- 
ment by Chiao et alO More recently, this method has 
been used to calculate a variety of quantities, which in- 
cludes the angular dependence of the specific heat_D|Lc?- 
wave superconductivity in a planar magnetic field.EjO 

The object of this paper is to study both the thermo- 
dynamics and the thermal conductivity tensor in a pla- 
nar magnetic field. We find the analytic expressions of 
the density of states and the thermal conductivity tensor 
in the presence of a magnetic field and impurity. Also, 
the Doppler shift is handled consistently with the layered 
structere common to both high-Tc cuprates and k-(ET)2 
saltsEJ This enables us to analyze these quantities both 
in the superclean limit (■£ ^ ^1) a.nd in the clean 
limit (-^ < ;£ < 1), while in [16-18] they are called the 
clean limit and the dirty limit. We think the word "dirty 
limit" should not be used to the system where ■£ ^ 1. 
Here Hc2 is the upper critical field, F the quasi-particle 
scattering rate, and A the order parameter. We find 
the angular dependence of both the specific heat and the 
superfluid density roughly by a factor 4 less than those 
predicted in [23]. Also, the longitudinal thermal conduc- 
tivity increases linearly with H in the superclean limit, 
while H \n{Hc2/ H) in the clean limit. Therefore, our 
result is similar to [18] in the clean limit, though it is dif- 
ferent from [18] in the superclean limit. So the thermal 
conductivity increases with H as in [18] but in contradic- 
tion to [13]. Further the thermal conductivity exhibits 
a similar ^-dependence as observed in [13], but again it 
has the opposite sign; Ky takes the minimum value for 
6 = zLj contrary to [13]. Perhaps this is due to the fact 
there is more quasi-particle for 9 = than for 9 = ±j. 
On the other hand, the transverse thermal conductivity 
is described by k± ~ —Hsm{29) in the both limits. This 
dependence appears to describe quite well the data re- 
ported in [11,12]. For simplicity, we limit our analysis 
to T <C Tc and H <C ffc2 ancjUjjb^ impurity scattering is 
treated in the unitarity limit. I 
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II. DENSITY OF STATES, SPECIFIC HEAT, AND 
SUPERFLUID DENSITY 

Following Barash et al,lii the impurity-renormalized 
quasiparticle energy Cj is given by 



tl) = cij + i 



and 



- (■ 



V • q 



V • qj 



A2 cos2(20) 



(1) 



(2) 



where |v 



I =^ vv' (sin^ X + sin^ {(/)- 9)), the 
Doppler-shifted energy due to the circulating supercur- 
rent around the vortex, r the distance from the vortex, 
V and v' are the Fermi velocity within the a-b plane and 
the one parallel to the c-axis, respectively, where the et. 
feet of the layered structure is explicitly considered .Eil 
Here x — cps, is the angle the quasi-particle mo- 
mentum in the a-b plane makes from the a-axis, 9 is 
the angle the magnetic field makes from the a-axis, and 
< . . . >= ^^^y2 J dxj d(f> . . . the angular average. We as- 
sume the unitarity limit of the impurity scattering. The 
low temperature limit of both the specific heat and the 
superfluid densiijj is easily obtained from the residual 
density of statesEII (i.e. the density of states on the Fermi 
surface or at w = 0). The quasi-particle density of states 
at w = in the presence of both impurities and a mag- 
netic field is given by 



N(uj = 0) 



Re g{uj = 0) 

7r\ 'V^TP 



)+xtan-i(^)) (3) 



where x = |v • q|/A and we used 



— I = T = 



(4) 



This set of equation is solved in the following: 
For Co « (x) « 1 (i.e. £ « « 1), 



N{lo = 0) 



N{H,i 



No ' 

-±^^(ln(l)-l), 
TT A (a;) X 



(5) 



and 



Co 



F 1 

A(^ 

H 



+ 



H., « ^ « 1)' 



For (x) < Co < 1 (i.e 
N{uj = 0),_ NiH,9)._ iVimp(O) 



1 A 



!^IIP^(1 + ±^1(2;2^) (6) 



and 



Co 



2 A 



/ln(4 



2 A 



where N^np (0) is the density of states in the H — P-iCase 
with the unitarity impurity scatterers and is givenc3 



iVo TT Go 



(7) 



We call the former the superclean limit while the latter 
the clean limit. These have been obtained essentially in 
[16] except for a few typos. We used the same angular 
and-thp spatial average as in Barash et al and Kiibert et 
al.E£lEZl Finally, the density of states in the both limits is 
the following: 

For the superclean limit (Cq <C (x) <C 1); 



N{H,9) Vvv'eH 



No 



where 



I{9) 



[ln( 



8A 



y/vv'eH 
(8) 



1(9) ^^J:±U sin' X + sin\±^ 



1, 



3±s 



3±s 



0.0285 cos(46') -I- 0.955 



(9) 



and 



Ji9) 



ii;±<^ln Jsin^x + sin2(±^ - 



^1 2±,s + V(3±g)(l±g) 

4 ^ 4 
~ -0.778 + 0.744 Max{ I cos 6l|, |sin6l|} (10) 

where s = sin(20) and E{k) is the complete elliptic in- 
tegral. Here, instead of average over 0, wej-average over 
(f> = ±j (i.e. over the two nodal directionsEj) 

The 6'-dependence of I {9} and J (9) are shown in Fig.l 
and Fig. 2. In Fig.l we compared also the angular de- 
pendent obtained in [23] where the layered structure is 
ignored. In the present calculation, the angular depen- 
dent term is about the 4% of the coefficient ~ Vh, while 
in [23] it is about 20 %. A recent specific data from sto- 
ichiometric YBa2Cu307.oo crystal^3 appear to be more 
consistent with the present analysis. Perhaps we have to 
point out the angular dependence of J{9) appeared since 
we cut off the the short-range logarithmic divergence at 
X = ^ where the uniform order parameter started to 
be greatly disturbed,c3 though the exact value ^ in the 
present case is not so important. 
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For the clean limit ((a-) < Cq < 1) 



N{H,e) ^ jVi.np(0) 
No No 

where 

m 



[1 + 7 



1 vv'eH , 



4 TA 



-(ln( 



vv'eH 



m)] 
(11) 



-S±^(sin^X + sin^(±^-^)) x 



In 



sm^ X - 



sin2(±- 



(2 ± s) ln( 



2±s + V(3±s)(l±s), 



4- V(3±s)(l±s) 
0.147- 0.082 cos(46i) 



(12) 



The function F{9) is shown in Fig. 3. 

The residual density of states is most readily accessible 
to the low temperature limit of the spin susceptibility as; 
in NMR, the specific heat and the superfluid density.EZi 
In particular the low temperature specific heat and the 
superfluid density are given by 



27r2 

^TNiH, 



and 



Ps{H,e) = l-N(H,9)/No 



(13) 



(14) 



Also, in the superclean limit and T/A > (x), we obtainEHl 

(15) 



= 18C(3) JtVo 



and 



p,(r)M0) = l-2(ln2) 



A 



(16) 



To conclude this section, we have considered the resid- 
ual density of states in a planar magnetic field both in 
the superclean and in the clean limit. In both limits the 
residual density of states exhibits the fourfold symme- 
try, though magnitude of thiSj-terms is roughly 4 times 
smaller than the earlier result E3 



III. THERMAL CONDUCTIVITY TENSOR 

Making use of Ambegaokar-Griffin formulalll the low 
temperature thermal conductivity tensor is given by 



ku/ko 



1 

TT / 2 

2 



"2-COs2( 
\2 , „„-2f 



■ X- - cos ^(20) \ 



IjCo+ixY +cos\2(j))\J 
+ ixY + cos2(2(/i) 



(17) 



Similarly 



1 / Cg + a;^-cos^(20) \ 
''(zinili) '^^ ^ |(Co-^^x)^-^co?(2^y 
2\ ^ ' ReV(<^o+ia;)2+cos2(20) 



(18) 

where kq — n\\ {H = 0) 

First, for the superclean hmit (Co ^ (a;) <C 1) Eq.(17) 
and Eq.(18) reduce to 



" TT 



and 



TT A2 



(19) 



(20) 



where k„ — , the thermal conductivity in the nor- 

dFto 
mal state, and 



1 



3 + /^) 



TT - V 2 

~ 0.29 sin(26') (21) 

The function L{9) is shown in Fig. 4 together with the 
approximate form. We note also that K|| has the same 
angular dependence as [N{H,d)]'^. Also k± is propor- 
tional to sin(20) in a good approximation (see Fig. 4 and 
Eq.(9)). 

On the other hand, in the clean limit((a::) ^ Co <C 1), 
we obtain 

K\\/Ko — i- + t; 



3 Ci 



= 1 



37r FA V V ttF^L ^V^^TIff 



m 



) 

(22) 



where F{9) is given in Eq.(12). Note that F{9) appears 
from the lower cut-off of r as explained after Eq.(lO). 
And 



k±/ko = 



1 vv'eH I ,^ /2A, r . . 2A , 

Svr FA V V ttF L \/vv'eH 



G{9) 



(23) 



where 



cw4[(2 + ,)m1±i±^S±311±3) 



(2- s) ln(- 



s+ V(3-s)(l-s). 



v/(3 + s)(l + s) + v/(3-s)(l-s) 
0.422 sin(26l) 



2 sin(26') 



(24) 
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The function G{9) is shown in Fig. 5 with the approximate 

form. Here, kq = — - — , n is the quasiparticle density 
3, Am 

and m the quasiparticle mass. This form of thermal con- 
ductiviUi. in the absence of magnetic field is derived first 
by Lee.Ea In general, however, A in kq depends on both 
on r and H. In the superclean limit, A(_ff) may be ap- 
proximately given hyS3 



A(H)/A„ 



for Co < (a;) < 1 



1 vv'eH 
6^2- 



In the clean limit ((x) ^ Cq <C 1), 

A(r)/Aoc.i-5l 



(25) 



(26) 



In both limits, the longitudinal conductivity increases 
with H, linearly in superclean limit and H\n{Hc2/H) 
in clean limit, respectively. (Our result agrees with [18] 
only in the clean limit, <C ■£ <C 1.) Further it ex- 
hibits the 9 dependence in both the superclean and the 
clean limit. In the superclean limit the 0-dependence 
comes from that of the density of states, while in the 
clean limit this arises from the short range cut-off we 
have introduced after Eq.(lO). Also this angular depen- 
dence is very similar to the one reported in [13], but of 
opposite sign. But perhaps of particular interest is the 
transverse thermal conductivity. The dominant terms in 
Eq.(20) and Eq.(23), exhibit sin(26') dependence, which is 
fully consistent with the early experiment. MO Also tlji 
sin(20) dependence is appreciable even when H ^ Hc2 



IV. CONCLUDING REMARKS 

We have extended earlier analysis of the thermody- 
namics and the transport properties in d-wave super- 
conductors in 2 directions. First, we take account of 
the layered structure of the underlying superconductors 
explicitly. Second, we focused on the angular depen- 
dence of thermal conductivity tensor, which exhibits clear 
signs of the nodal structures in d-wave superconductors. 
Indeed the diagonal thermal conductivity exhibits the 
fourfold symmetry as observed in [13], but of opposite 
sign. In particular, the present result describes the trans- 
verse thermal coBd|ttctivity (or Righi-Leduc effect) ob- 
served in YBCO.EBEj Indeed, we have found recently a 
similar transverse thermal conductivity in p-wave super- 
conductocs as Sr2Ru04 and /-wave superconductors as 
in UPts.Ej Therefore, this sin(20) dependence in k± is 
rather common to most of unconventional superconduc- 
tors. 

In summary, the exploration of the nodal structure in 
unconventional superconductors will provide useful in- 
sight in the quasi-particle spectrum in the vortex state. 
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FIG. 1. I (9) and the approximate form are shown as func- 
tion of 6. Ii{6) is the angular dependence in [23] where the 
layer structure is ignored. 



FIG. 2. J{6) is shown as a function of 6. It has cusps at 
7r/4 and 37r/4 together with the approximate form. 

FIG. 3. F{9) is shown as a function of 9 together with the 
approximate form. This angular dependence appears in the 
specific heat Cv, the superfiuid density ps, and thermal con- 
ductivity tensor in clean limit(-^-j <^ ^ 1). 

FIG. 4. L[9) is shown as a function of 9 to together with 
the approximate form. 

FIG. 5. G{9) is shown as a function of 9 together with the 
approximate form. 
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